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Abstract 

The theoretical investigations on electrorheological (ER) fluids are usuahy 
concentrated on monodisperse systems. Real ER fluids must be polydisperse 
in nature, i.e., the suspended particles can have various sizes and/or different 
dielectric constants. An initial approach for these studies would be the point- 
dipole (PD) approximation, which is known to err considerably when the 
particles approach and finally touch due to multipolar interactions. In a 
recent work, we proposed a dipole-induced-dipole (DID) model for computer 
simulation of ER fluids, which was shown to be both more accurate than 
the PD model and easy to use. The DID model was applied to simulate 
the athermal aggregation of particles in ER fluids and the aggregation time 
was found to be signiflcantly reduced as compared to the PD model. In this 
work, we will report results for the case when the dielectric contrasts of some 
particles can be negative. In which case, the direction of the force is reversed. 
Moreover, the inclusion of DID force further complicates the results because 
the symmetry between positive and negative contrasts will be broken by the 
presence of dipole-induced interactions. 
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I. INTRODUCTION 



Many theoretical investigations on electrorheological (ER) fluids are usually concentrated 
on monodisperse systems in which all the suspended particles are of the same size and di- 
electric constant. In reality ER fluids must be polydisperse in nature, i.e., the particles can 
have various sizes and/or different dielectric permittivities For instance, the particle 
size has a significant impact on the yield stress [Q], as well on the rheology An initial 
approach for these studies would be the point-dipole (PD) approximation As many-body 
and multipolar interactions between the particles have been neglected, the PD approxima- 
tion is known to err considerably when the particles approach and finally touch. The PD 
approximation becomes even worse when the dielectric contrasts between the suspended 
particles and the host medium become large. To circumvent the problem, we recently em- 
ployed the multiple image method to compute the interparticle force for a polydisperse ER 
fluid 0. From the results, we proposed a dipole-induced-dipole (DID) model for computer 
simulations of ER fluids; the DID model yields very good agreements with the multiple 
image results for a wide range of dielectric contrasts and polydispersity [^]. The DID model 
has recently been employed to simulate the athermal aggregation of particles in ER fluids 
in which the particles are of different permittivities 0. Moreover, the dielectric contrasts 
between the particles and the host fluids were all positive. The aggregation time was found 
to be significantly reduced, both in uniaxial and rotating electric fields j^]. In this work, 
we will report results for the case when the dielectric contrasts of some particles can be 
negative. In which case, the direction of the force should be reversed [^. The inclusion of 
the DID force further complicates the results because the symmetry between positive and 
negative contrasts will be broken by the presence of dipole-induced interactions. We found 
that the aggregation time can be much increased as compared to the PD model. 

In the next section, we review the multiple image method and establish the DID model. 
In section III, we apply the DID model to the computer simulation of ER fluids in a uniaxial 
field. In section IV, we extend the simulation to athermal aggregation in rotating fields. 
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Discussion on our results will be given. 



II. MULTIPLE IMAGE METHOD 



In this section, we generalize the multiple image method |^ to handle both positive and 
negative dielectric contrasts. Consider a pair of dielectric spheres, of radii a and 6, dielectric 
constants ei and e'^ respectively, separated by a distance r. The spheres are embedded in 
a host medium of a dielectric constant €2- Upon the application of an electric field Eq, the 
induced-dipole moment inside the spheres are, respectively, given by (SI units) 
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where the dipole factors P, /5' are defined as 

ei - £2 



ei - £2 



ei + 2e2 e'i + 2e2 

In the point-dipole (PD) model, the force between two particles is given by 
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where F = 12'n'eoe2P (3' a^b^ , r and 6 are unit vectors. Many previous studies were concen- 
trated on the case in which both /3, /?' > 0. When the /3s adopt opposite signs, the direction 
of the PD force will be reversed. In a recent paper, we derived a correction to the PD force 
from the multiple image method H. The total dipole moment inside sphere a is 
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where the subscripts T(L) denote a transverse (longitudinal) field, i.e., the applied field is 
perpendicular (parallel) to the line joining the centers of the spheres. Similar expressions 
for the total dipole moment inside sphere b can be obtained by interchanging a and 6, as 
well as (3 and (3' . The parameter a satisfies 
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The forces between the spheres is given by 0] 
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It should be noted that the multiple image results can be used to compare among the 
various models according to how many terms are retained in the multiple image expressions: 
(a) PD model: n = 1 term only, (b) DID model: n = 1 to n = 2 terms only, and (c) 
multiple-induced-dipole (MID) model: n = 1 to n = oo terms. 

For convenience, we define the reduced separation a = r/{a + b). Here we set a = b. We 
consider two cases: (a) /5 > 0, /3' < or ei > 62 > €[ and (b) /5 = < or £2 > Ci = e[. The 
interparticle forces in the longitudinal and transverse cases are plotted in Fig.|l| with different 
values of the dielectric contrasts. For the ) depicted in Fig.|l|(a), the magnitude of 

MID force falls between that of the PD and DID in both the longitudinal and transverse 
cases. At low contrasts, the DID results almost coincide with the MID results. While at 
high contrast, the DID model exhibits significant deviation from MID when a < 1.1. For the 
case (b) as depicted in Fig.^(b), the DID results almost coincide with the MID results at low 
contrast, while deviate significantly at high contrast. The forces are found to be qualitatively 
different from the case of /3, f3' > 0. Thus, by including the multiple image contributions, it 
is observed that the case of (3, (3' < significantly different from the presumably symmetric 
case of (3,13' > 0. The symmetry has been broken due to the presence of DID forces. 

III. ATHERMAL AGGREGATION IN THE UNIAXIAL FIELD 

The multiple image expressions [Eqs.(0) and (1)] allow us to calculate the correction 
factor defined as the ratio between the DID and the PD forces [^: 
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where = 3paoPbo/4vreoe2r^, F^J'^ = -6paoP6o/47reoe2r^, and F^^ = -3paoPw/4vreoe2r^ 
are the point-dipole forces for the transverse, longitudinal, and F cases, respectively. If we 
denote the ratios in Eqs.(|])-(|l^) by Ki\,K± and Kr respectively, the force between two 
particles is modified to 

F 
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In what follows, we consider two spheres of equal radius a and various dipole factors f3 
and (3' respectively. These spheres are initially at rest and at a separation d^. For aggregation 
induced by a uniaxial field, we consider two cases. For case (a) /5 > 0, < or ei > e2 > e[, 
the electric field is perpendicular to the line joining the centers of the spheres. The equation 
of motion is given by 

|^Fx(2.), (11) 

while for case (h) P = P' < or ei = e'l < 62, the electric field is parallel to the line joining 
the centers of the spheres. The equation of motion becomes 

1-^,(2.), (12) 

where z is the displacement of one sphere from the center of mass. The separation between 
the two spheres is therefore d = 2z and the initial condition is d = do at t = 0. Eqs.(|ll]) 
and (|12D are dimensionless equation. Following Klingenberg and with slight modification, 
we choose the following natural scaling units to define the dimensionless variables 0: 
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zo = 2a, to = UTTricO^/Fo, Fq = -7^eoe2Fo«^ (13) 

where Eq is the field strength, m the mass, and rjc the coefficient of viscosity. We have 
followed Klingenberg to ignore the inertial effect and thermal motion of the particles [Q. 
The initial separation d^ is related to the volume fraction (p by 



In the PD approximation, Eq.(|rT]) admits an analytic solution 
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while Eq.(]T^) gives 
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It is obvious the above equations impose a condition for aggregation in uniaxial field, namely, 
the transverse case requires that f3 and (3' are of different signs while the longitudinal case 
requires that they are of the same sign. 

For the DID model, we integrate the equations of motion by the fourth-order Runge- 
Kutta algorithm, with the time step 6t = 0.0001 for both small and large volume fractions. 
In Fig.^ we plot the reduced separation a against dimensionless time t/to- The results reveal 
that in general the DID results deviate slightly from the PD results at low volume fractions, 
i.e. at a large initial separation, while the deviation becomes large at high volume fractions. 
The deviations are more pronounced at high contrasts, attributed to a large attractive force, 
resulting in a smaller aggregation time. In both cases (a) and (b), the DID aggregation time 
is generally larger than that of PD. While in the case (3,13' > 0, the situation is reversed. 
This feature can be understood from the force magnitudes of the two models (see Fig.|l]). 

Fig.0 shows the ratio of the aggregation time for the PD model to the DID model against 
reduced initial separation do /2a. The aggregation time of the DID model is generally larger 
than that of the PD model especially when the initial separation is small. The correction 
factor is more pronounced at high contrast. In the case /3 = < or 62 > ei = e'^, 
we observe a non-monotonic dependence, the correction factor increases significantly when 
the spheres are close and decrease again when d/2a = 1.13 or less at high contrast. This 
behavior can be understood from the interparticle force (Fig.|I](b), the three panels on the 
right). At low contrast, the difference in magnitude between the PD and DID force increases 



monotonically as the separation decreases. However at high contrast the difference increases 
first then decreases again as the separation decreases. 



IV. ATHERMAL AGGREGATION IN THE ROTATING FIELD 

Perhaps it is more interesting to consider aggregation in a rotating electric field [Q. 
Consider a rotating field applied in the x-y plane, = Eocosut, Ey = Eosmurt. The 
dimensionless equation of motion for the two spheres becomes 

— = F|| (r) cos^ cut + F±{r) sin^ cut, -j- = -Fr{r) sin 2ujt (17) 

where (x, y) are the displacement of one sphere from their center of mass and r = 2A/x^~+~y^ 
is the separation between two particles. In case of a large u, we may safely neglect the 
y direction of the motion. For the PD approximation, the dimensionless forces are Fy = 
—2[3(3'/r^ and Fj_ = [3[3' /r^, respectively, which yields the analytic solution 
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The separation between two spheres is d = 2x. For the DID model, we integrate the equation 
of motion by the fourth-order Runge-Kutta algorithm, with a time step 6t = l/400u;. 

Note that the two particles cannot aggregate if /? > 0, /3' < or ei > £2 > e[ in the 
rotating field because the repulsive force induced by longitudinal field F^ is always larger 
than the attractive force induced by transverse field Ft, while the particles spend equal 
times in both fields on the average. The displacement-time graph of the aggregations with 
62 > ei = e'l in the two models are plotted in Fig.^. The shape of the displacement- 
time graph of DID model is different from that of the PD model, especially when the two 
spheres are close, resulting in a large deviation between the DID and PD models in the 
aggregation time. This behavior can again be understood by the interparticle force between 
the two spheres. We can see from Fig.|l|(b) that the DID longitudinal attractive force and its 
transverse repulsive force is about the same when the separation is small. Since the spheres 
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spend equal times in the longitudinal and transverse field and they change the direction of 
motion according to the rotating field, the displacement of the spheres can be estimated by 
the difference between the magnitude of the longitudinal attractive force and the transverse 
repulsive force, as there would be no displacement if the magnitude of the two forces are 
the same. Thus when the two spheres are near, their velocity become smaller and have 
different shape in the displacement-time graph. On the other hand, the PD attractive force 
is sufficiently larger than its repulsive force at any separation. 

We further calculated the correction factor of the aggregation time for the PD model with 
respect to the DID model in the same way as in the uniaxial field case. We plot the correction 
factor against the reduced initial separation d^jla with different parameters in Fig.^. The 
aggregation time is significantly increased when the mutual polarization of the two spheres 
is taken into account. The increase becomes even more pronounced and oscillating when the 
initial separation is small. The oscillation is due to the sensitive dependence on the initial 
polarization of dipoles when the spheres are close. 

DISCUSSION AND CONCLUSION 

Here a few comments on our results are in order. In this work, we have concentrated on 
the multipolar interactions between touching particles, which we believe are more important 
than the many-body or local- field effects In the latter approaches, the particles in 

ER fluids are still treated as point dipoles, while their dipole moments are determined 
by adding the local-field corrections. In the DID model, the additional terms arise from 
multipole interactions, rather than from local-field corrections. 

We have studied the aggregation time for two spherical particles. We should also examine 
the morphology of aggregation in polydisperse ER fiuids, due to the dipole-induced forces. 
To this end, it may be difficult to obtain the ground state lattice in polydisperse ER fiuids 
by using dynamic simulations alone, since both body-centred tetragonal and simple cubic or 
even honeycomb stacking lattices coexist in such systems according to Kaski and coworkers 
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TT| . Moreover, many external factors such as the numbers of particles, the size of the cell, 
etc. will affect the final results. We are currently trying to incorporate the DID terms in the 
interaction energy. In this connection, we can also examine the recently proposed structural 
transformation by applying the uniaxial and rotating fields simultaneously to an ER fluid 

In summary, we have used the DID model to deal with computer simulations of polydis- 
perse ER fluids for the case when the dielectric contrasts of some particles can be negative. 
We studied athermal aggregation of two spherical particles both in uniaxial and rotating 
electric fields. We showed that an inclusion of the DID force breaks the symmetry between 
positive and negative contrasts. As a result, the aggregation time can be much increased as 
compared to the PD model. 
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FIGURES 

FIG. 1. The interparticle force of the PD and DID models plotted against the reduced sepa- 
ration a between two spherical particles for several dielectric contrasts, ei/e2 for one particle and 
e'i/e2 for the other, both in the transverse (T) and longitudinal (L) electric fields: (a) /? > 0, < 
or ei > 62 > e'l and (b) /3 = < or €2 > ei = e^. 

FIG. 2. The displacement-time graph for athermal aggregation of two spherical particles in a 
uniaxial field for various dielectric contrasts: (a) /? > 0, /9' < or ei > 62 > e'^, (b) €2 > ei = e'l in 
low volume fractions, and (c) €2 > ei = e'-^ in high volume fractions. 

FIG. 3. The correction factor of the aggregation time of two spherical particles in a uniaxial 
field plotted against the initial separation. The upper panel corresponds to the transverse field 
case while the lower pannel corresponds to the longitudinal field case. 

FIG. 4. The displacement-time graph for athermal aggregation of two spherical particles in a 
rotating field for various dielectric contrasts, rotating field frequencies for (a) low volume fractions 
and (b) high volume fractions. 

FIG. 5. The correction factor of the aggregation time of two spherical particles in a rotating 
field plotted against the initial separation for two different rotating field frequencies w, and tdid 
and TpD are the aggregation times in the DID and PD models respectively. 
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Fig. 1(a) /Wong, Sun and Yu 
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Fig. 2 (a) /Wong, Sun and Yu 
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Fig.3/Wong, Sun and Yu 
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Fig.4(a)/Wong, Sun and Yu 
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